Quantifying structural damage from self-irradiation in a plutonium superconductor 
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The 18.5 K superconductor PuCoGas has many unusual properties, including those due to dam- 
age induced by self-irradiation. The superconducting transition temperature decreases sharply with 
time, suggesting a radiation-induced Frenkel defect concentration much larger than predicted by 
current radiation damage theories. Extended x-ray absorption fine-structure measurements demon- 
strate that while the local crystal structure in fresh material is well ordered, aged material is dis- 
ordered much more strongly than expected from simple defects, consistent with strong disorder 
throughout the damage cascade region. These data highlight the potential impact of local lattice 
distortions relative to defects on the properties of irradiated materials and underscore the need for 
more atomic-resolution structural comparisons between radiation damage experiments and theory. 

PACS numbers: 71.27. +a, 74.70.Tx, 61.80.-x, 61.10.Ht 
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I. INTRODUCTION 

Plutonium is arguably the most complex and least un- 
derstood of all elements, ultimately due to the propensity 
of its 5 /-electrons to simultaneously reside in bonding 
and non-bonding electronic states. Theoretical models 
of S- and cv-Pu, for instance, show promise for explaining 
this complex behavioral 3 - but assume a homogeneous 
crystalline structure despite unavoidable self-irradiation 
damage, which can significantly alter magnetic and elec- 
tronic properties. Superconductors provide a path for 
elucidating radiation damage effects^ since their proper- 
ties are especially vulnerable to atomic-level structural 
disorder. Here, we report extended x-ray absorption 
fine-structure (EXAFS) measurements on the PuCoGa 5 
superconductor— that demonstrate the local structure of 
aged material is damaged at least an order of magnitude 
faster than theoretical predictions focusing on Frenkel 
defects indicated These results explain the sharp reduc- 
tion of the superconducting critical temperature, T c , with 
time and underscore the need for improved radiation 
damage models relevant not only for understanding plu- 
tonium superconductors, but also for Pu metal and other 
radioactive materials. 

In addition to changing fundamental properties in ele- 
mental Pu and PuCoGa,5, radiation damage affects many 
aspects of science and industry, for example, in semi- 
conductors, nuclear power generation and its associated 
waste disposal, and the aging nuclear stockpile. Con- 
sequently, radiation damage has been studied for over 
a century, and intensely since World War II. Although 
experimental measures have provided ample verification 
of long-range (>10 nm) structural effects due to radia- 
tion damage, atomic-resolution descriptions of the dam- 
age induced in a crystal have relied almost exclusively 
on theoretical calculations. This reliance has been pre- 
dictably dangerous. Recenti^ nuclear magnetic reso- 



nance (NMR) experiments reveal that radiation damage 
accumulates about 5 times faster in zircons and other 
ceramics than calculations indicate, calling into question 
both the theoretical calculations and the viability of nu- 
clear waste containment schemes. Experimental stud- 
ies have not generally kept pace with theoretical treat- 
ments of fully relaxed damage cascade structures at the 
atomic level. Damage cascade calculations have not been 
quantitatively compared to any atomic-resolution exper- 
iments until very recently^ and never to our knowledge 
in intcrmctallics 9 or technologically important materials 
such as <5-Pu. Measurements of such structural changes 
are vital to verify and improve these theories, and for 
comparison to more complex theories of dislocation loops, 
He bubble formation, and volume expansion. 

Because large regions of a given sample remain in an 
undamaged state while damaged regions no longer have 
translational symmetry, traditional scattering techniques 
such as Rietveld analysis of powder diffraction data have 
not generated a quantitative measure of the damage frac- 
tion of a material. Local probes that treat ordered and 
disordered regions on an equal footing should be able to 
provide a more detailed, atomic-level description of the 
damage. NMR is one such technique, although it does 
not easily lend itself to a direct structural interpretation. 
The EXAFS technique employed here provides another 
local probe and has the advantage that it gives radial 
pair-distance distribution information around a specific 
atomic species, since it relies on the backscattering of a 
photoelectron from a core excitation. 

To investigate structural radiation damage effects at 
the atomic level in plutonium, we employ superconduct- 
ing PuCoGas (Fig. [T|). The unusual properties of this 
18.5 K superconductor, i.e., a nearly ten- fold higher T c 
than any other heavy-fermion related intermetallic and 
an electronic structure that bears a strong resemblance to 
5-Pu, have been well characterized ^ 10 ' 11 Superconduct- 
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FIG. 1: (Color online) The superconducting transition tem- 
perature (T c ) as measured by magnetometry as a function of 
the a-decays per atom for this study (samples A and B) and 
in the literature (1 yr=4.9xl0~ 6 a-decays per atom for these 
samples). A linear decay of T c ( — ) starting from about 19.0 
K and decaying at a rate of about 5.4 x 10 5 K per a-decay per 
atom is also shown for reference. The inset shows the crystal 
structure of PuCoGas .— 



ing properties that change with time due to plutonium 
a-decay include a decrease in T c by ^0.2 K per month 
(Fig. [T|). the upper critical field H C 2 ~ 70 T, and large 
critical current density J c > 10 4 A/cm 2 for T > 0.9T C ^ 
likely caused by damage-induced scattering and pinning 
centers. In addition, the impurity scattering rate inferred 
from NMR experiments follows the observed change in 
T c ^ 

Although a modern theoretical treatment of radiation 
damage in PuCoGas is currently not available, we utilize 
calculations on the structurally similar i5-Pu system for 
comparison^ In these models, the a-particle generated 
by the decay of a 239 Pu nucleus has about 5 MeV of en- 
ergy and ballistically generates ~300 Frenkel defect pairs 
over a distance of nearly a micron^ Most of the dam- 
age, however, is done by the recoiling 235 U nucleus with 
86 keV, which produces ~2300 Frenkel pairs. A typical 
generated damage cascade extends over nearly 10 nm, 
with a defect volume fraction of about 3%. It is impor- 
tant to note that only the effect of Frenkel-type defects 
are considered in this view of radiation damage; how- 
ever, these models form the basis of most modern theo- 
ries such as molecular dynamics theories (for example, see 
Ref. [l~2l ). In fact, lattice relaxation can occur during the 
intermediate time-scales after these displacement events, 
and have been explored in various metallic systems with 
molecular dynamics and kinetic Monte Carlo techniques. 
These latter calculations show that the effective number 



of defects is reduced by as much as an order of mag- 
nitude within only a few picoseconds due to additional 
defect migration! 1 ^ These values should not change sub- 
stantially in PuCoGa.5, and in fact a rough TRIM code 1 ^ 
calculation (which does not include the longer-time scale 
defect migration) using default values generates similar 
damage rates and ranges to the models above. We there- 
fore expect far fewer than the initial A/d ^2600 pair de- 
fects to survive the defect migration per a-decay. Using 
an a-decay rate X a w 3.4 x 10~ 5 per Pu per year from the 
samples discussed below, the expected upper-limit dam- 
age fraction / tot = 2Ni)X a /7 f=a 2.5% after one year, with 
l/fth f the atoms being Pu. T c is reduced by 50% after 
about 3 years (Fig. [T]). These damage estimates therefore 
indicate that, including the expected additional defect 
migration, f tot (T c /T c0 = 50%) < 7.5%. This estimate 
is inconsistent with what one expects in a short coher- 
ence length superconductor: Within the Bogoliubov-de 
Gennes (BdG) formalism with strong scattering and a 
coherence length £0 w 2.1 nm£), the lower limit (strong 
scattering limit) damage fraction is fBdG(T c /T c o = 
50%) «15% J£ This value is consistent with recent studies 
of Cei-zLazCoIng 1 ^ and CeCoIns^SnaJi where the sub- 
stitutions are mostly within the superconducting planes, 
and therefore likely produce strong scattering. The crit- 
ical damage fraction should be higher when defects are 
randomly distributed, as with radiation damage. 

This paper continues with a description of the sample 
characterization and other experimental methods (Sec. 
HI)) , details of the data analysis and results (Sec. IIIip . a 
discussion of the implications of these results (Sec. IIV|) 
and a concluding summary (Sec. |V| . 



II. EXPERIMENTAL METHODS 

Two PuCoGas samples were synthesized from Ga 
flux.— At the time of the most recent x-ray measure- 
ments, one was about 3 years old (sample A), and the 
other was about 1 year old (sample B). Superconduct- 
ing critical temperatures were measured in a Quantum 
Design Magnetic Properties Measurement System as the 
point at which diamagnetism was observed in 10 Oe (Fig. 
[TJ. The isotopic content is the same for both samples 
with the main radioactivity coming from 93.93% 239 Pu, 
5.85% 240 Pu, and 0.12% 241 Pu. The accumulated dose 
has been shown to be a reasonable indicator of how sam- 
ple properties change with time^ and we report the sam- 
ple age in units of a-decays per atom (1 yr=4.9xl0~ 6 a- 
decays per atom in the formula unit). We do not use the 
more common "displacements per atom" (dpa) unit, be- 
cause it generally assumes the number of displaced atoms 
(Frenkel pairs) is known from theoretical models, and we 
show below that this assumption may not be correct. 

The material was triply contained for EXAFS ex- 
periments using epoxy- and indium-scaled kapton win- 
dows, and placed into a LHe flow cryostat at T <30 K. 
Sample A was initially measured (measurements with 
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FIG. 2: (Color online) EXAFS data k 3 X (k) vs. k from the 
(a) Pu Lin edge, (b) the Co K edge, and (c) the Ga K edge, 
for a fresh ( — , 0.2 x 10 -6 a-decays per atom) and an aged 
( — , 14.5 xlO -6 a-decays per atom) sample. 



< 4 x 1CP 6 a-decays per atom) as a single crystal in 
fluorescence mode using 30-element Ge detectors, with 
the data corrected for dead time and self absorption^ 
Otherwise the samples were ground for the EXAFS ex- 
periments and passed through a 32 /im sieve, with about 
8 mg of this powder mixed with dried boron nitride and 
packed into an aluminum frame. Transmission and fluo- 
rescence data agree quantitatively. This sample mass re- 
sulted in a change in absorption across the Pu Ini edge 
Afi a of ~0.5 absorption lengths, whereas A/i a ~0.3 and 
A^i a ~2.2 absorption lengths across the Co and Ga K 
edges, respectively. EXAFS spectra were collected at the 
Stanford Synchrotron Radiation Laboratory on beam- 
lines 10-2 and 11-2 over a three year period, at the Pu 
Lui, Co K, and Ga K edges, generally using a half-tuned, 
double crystal Si(220) monochromator. The monochro- 
mator resolution was adjusted such that it was well be- 
low the core-hole lifetime at a given edge. The data were 
analyzed using standard procedures 3i In particular, the 
embedded atom absorption /xo was determined using a 
cubic spline with between 4 and 6 knots over the data 
range, which was typically about 1 kcV above the ab- 
sorption threshold, Eq, as determined by the position of 
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FIG. 3: (Color online) Fourier transform (FT) of the k 3 \{k) 
data of the fresh sample ( — ) in Fig. [3] together with a fit 
to these data (---). The fit quality is such that the fit is 
difficult to distinguish from the data. The outer envelope is 
± the transform amplitude and the inner modulating line is 
the real part of the complex transform. The Pu edge data 
(a) are transformed between 2.5-16.0 A -1 , Gaussian broad- 
ened by 0.3 A -1 , and are fit between 2.0 and 5.0 A. The Co 
edge data (a) are transformed between 2.5-11.0 A" 1 , Gaus- 
sian broadened by 0.3 A -1 , and are fit between 2.0 and 5.5 A. 
The Ga edge data (a) are transformed between 2.5-15.0 A -1 , 
Gaussian broadened by 0.3 A -1 , and are fit between 1.6 and 
5.0 A. 

the half-height of the absorption change at the edge. The 
data were fit in r-space using the RSXAP packaged with 
theoretical scattering functions generated by FEFF7i^ 



III. ANALYSIS AND RESULTS 

Figure [5] shows an example of the normalized oscil- 
lations in the absorption above each measured edge as 
a function of k, the photoelectron wave vector, for the 
fresh (0.2 xlO -6 a-decays per atom) and 3 year old 
(14.5 xlO -6 a-decays per atom) samples. A Fourier 
transform (FT) of such data (Fig. [3]) produces peaks 
in the amplitude as a function of the distance r from the 
absorbing atomic species corresponding to neighboring 
atoms. For instance, the dominant peak in the Pu edge 
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TABLE I: EXAFS fit results for the Pu Lni and the Co K edges on a fresh (0.2 x 10 -6 a-decays per atom) sample of PuCoGa 5 . 
Fit and transform ranges are listed in Fig. [3] All single-scattering peaks within the fit range are included in these fits. Multiple 
scattering was only included in the Ga edge fits to avoid errors originating from peak overlap, but the results from these 
scattering paths are in themselves unreliable and so are not reported. Coordination numbers, N, are held fixed to the nominal 
structure. Sg, AE , and R(%) are 0.89(5), -4.7(3) eV and 3.6% for the Pu edge, 0.85(5), 5(1) eV, and 6.1% for the Co edge, 
and 0.85(3), -1.7(4), and 6.9% for the Ga edge, respectively. The number of free parameters in the fits are 10, 14, and 15 for the 
Pu, Co , and Ga edges, respectively, and are far below the number of independent data points as given by Stern's rule.— The 
quoted error on each quantity is the greater of that obtained by comparisons to standard materials,— repeated measurements, 
and a Monte-Carlo method.— The nature of the Ga local environment required constraints on several parameters to obtain 
meaningful fits, as indicated. Constraints on the pair distances in the Ga edge fits assume a tetragonal 115 structure such that 
Ga(l) sits at the center of the Pu ab face and a plane of Co atoms splits a plane of Ga(2). 
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data (Fig. [3^) corresponds to the 12 Pu-Ga neighbors at 
^3.0 A. Note that the atomic scattering functions gen- 
erate complicated lineshapes, causing shifts in the peak 
positions from the actual structure that are well repro- 
duced by calculations using FEFF7.— Information about 
the local structure is therefore obtained by fitting these 
calculated scattering functions to these data. The results 
of fits to data from fresh material that has not under- 
gone a significant amount of a-decay are reported in Fig. 
[3] and Table Q] The fit quality is excellent, and all the 
measured pair distances agree well with diffraction re- 
sults. We note that the local environment around Ga is 
more complicated than around Pu or Co, and therefore 
constraints were necessary to reduce the number of fit 
parameters while still obtaining high quality fits. Unfor- 
tunately, constraints add an unknown amount of system- 
atic error that is not reflected in the estimated errors, and 
we ascribe discrepancies between diffraction and EXAFS 
results from the Ga -ftT-edge data to this source. Where 



the data are of sufficient quality, the temperature de- 
pendence of the mean-squared displacements of the pair 
distances, a s, were obtained and arc well described by a 
correlatcd-Dcbye models with reasonable values of the 
correlatcd-Dcbye temperatures, c d's, (for comparison 
to <5-Pu, see Ref. I2H ) , and no evidence of static disorder 
from the fitted offsets, of tat 's. 

Significant radiation damage effects are readily appar- 
ent in the raw data with a marked decrease in the overall 
amplitude of the spectra as samples are aged up to 3 years 
(Fig. d]). Preliminary fits to the data from the aged sam- 
ples showed that the decrease in amplitude in the data is 
due both to a decrease in the overall scale factor Sq and 
an increase in each atom pair's a 2 with age. This situa- 
tion indicates that there are at least three distinct regions 
within the aged samples: virtually undamaged, strongly 
damaged where the distance widths a 2 are large enough 
that the local structure no longer contributes to the EX- 
AFS amplitude, and mildly damaged where the c^'s only 
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FIG. 4: (Color online) FT of the k 3 x(k) EXAFS data from 
the (a) Pu Liu, (b) Co K, and the (c) Ga K edges are shown 
for sample ages about one year apart. Samples have accumu- 
lated doses of 0.2xl0 -6 a-decays per atom ( — ), 5.4xl0 -6 a- 
decays per atom ( — ), 10.2 xlO -6 a-decays per atom ( — ), 
and 14.5 xlO -6 a-decays per atom ( ). Transform ranges 
are between 2.5-16.0 A" 1 , 2.5-10.0 A" 1 , and 2.5-14.5 A" 1 , 
respectively, all Gaussian broadened by 0.3 A -1 . 



allow for a weak contribution. The latter regions may 
exist, for instance, on the edges of strongly damaged re- 
gions. We therefore describe / to t as due to the sum of 
the strongly damaged fraction / s and the mildly dam- 
aged fraction f m . Since EXAFS amplitudes A ~ 1/a, 
as long as <r m is large enough, / tot » 1 - S 2 '(t)/S 2 (0), 
where S 2 '(t) ^ s obtained from fits where a 2 (t) are fixed 
at c 2 (0). That is, the total damage fraction in aged sam- 
ples can be estimated by fixing most fitting parameters 
to those obtained from fresh-sample fits, and estimating 
/tot by the change in amplitude as given by S 2 . These 
damage fractions are shown in Fig. and demonstrate 
strong damage production at apparently different rates 
for each atomic species. 

Using this estimate of /tot, fa is approximated by per- 
forming a fit where the a 2 (t)'s are no longer constrained, 
in which case / s « (S 2 (t) - S 21 '(*)) / S 2 (0) . Correlations 
between the Sq and a 2 parameters are more difficult to 
control in such a procedure, but fits to Pu edge data in- 
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FIG. 5: (Color online) The total fraction of displaced atoms as 
a function of the number of a-decays per atom, as seen from 
each atomic species. Error bars are based on reproducibility. 
Absolute errors are about ±5%. The line shows percolation 
behavior based on the first two aged Pu data points for sample 
A, and the hashed area gives the theoretical range of damage 
up to the amount expected in the absence of defect migration. 



dicate that between 60-80% of the damage is due to the 
strongly damaged regions and tr m is enhanced by ~ 0.0 10- 
0.015 A 2 over the well-ordered regions for all of the aged 
samples measured. Similar results are obtained from the 
Co and Ga edge data. 



IV. DISCUSSION 

The accumulated damage after one year in sample A 
is measured as between 20-30% from the point of view of 
each of the constituent atoms in PuCoGa,5. This value is 
an order of magnitude higher than theoretical estimates 
that only account for defects, or closer to two orders of 
magnitude if one allows for defect migration after defect 
formation^ It is important to note at this point that the 
damage enhancement, while possibly indicating more de- 
fects than expected, may also indicate significantly more 
defect-induced lattice strain and distortions. In fact, this 
amount of damage is so large that, assuming a ~10 nm 
damage cascade as expected for (5-Pu, every atom in the 
cascade is displaced from its equilibrium position. Mod- 
ern molecular dynamics calculations should, in principle, 
generate significant lattice distortions around defects in- 
duced by radiation damage, but this effect has not, to 
our knowledge, been reported in the literature for inter- 
metallics. In addition, the observed damage does not 
proceed at as fast a rate as extrapolated from lower ac- 
cumulated doses. This disagreement is very likely due to 
self-annealing 2 ^ caused by the room-temperature storage 
of the samples between measurements. Self-annealing is 
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also likely the cause of the differences between the sam- 
ples (sample B appears to be damaging at a somewhat 
slower rate), as the exact history of the storage condi- 
tions for each sample then becomes important. In addi- 
tion, differences between the constituent atomic species 
indicate different defect production or migration rates. 
These issues should be studied in other materials to fur- 
ther explore the role of self-annealing and atomic-species 
effects. 

Another way of describing radiation damage is to con- 
sider that a decay event generates enough heat within 
the damage cascade that the material locally melts and 
then rapidly resolidifies, thus quenching disorder into the 
cascade region from the high temperature state and po- 
tentially also creating a distribution of competing struc- 
tural phases. This picture is similar to the "thermal 
spike" model as originally proposed by Seitz,— although 
it proved to be calculationally intractable^ Applying 
molecular dynamics after (or during) Frenkel defect pro- 
duction seems to combine the relevant aspects of these 
two schemes, and could, in principle, generate the sort of 
distortions measured here in PuCoGas. Simulated EX- 
AFS data could be calculated from the results of future 
molecular dynamics calculations on PuCoGas for a direct 
comparison of the efficacy of such models. 

These data largely explain the fast reduction of T c in 
PuCoGa 5 , with T c /T c0 w 50% when f tot w 40%. This 
value is now greater than the lower-limit damage frac- 
tion of /BdG(T , c/T c o = 50%) «15%r^ as expected since 
not all damaged regions will create strong scattering. 
In addition, it is likely that defects outside the super- 
conducting planes scatter more weakly than defects on 
the in-plane sitesi 16 ! 17 Decreasing the effective scattering 
strength would increase the necessary impurity densities 
in the theory. 

Also shown in Fig. [5] is the prediction of a cubic per- 
colation model. The time axis is chosen such that the 
model agrees with the first aged Pu edge data points 
for sample A. Because of the strong damage that occurs 
within a cascade, the superconducting fraction is likely 
more closely related to the fraction of the material that 
exists within the volume between the edges of damage 
cascades, which is better described by the percolation 
model depicted in Fig. According to this extrapola- 
tion one might expect superconductivity to cease between 
3.5 and 4 years for these samples (/tot ~ 0.8), roughly 
consistent with the data in Fig. [TJ although no samples 
have yet been observed to become non-superconducting 
at this time. This simplification, of course, doesn't ac- 
count for any proximity effects, which should increase 
this time period, or any increased impurity scattering, 
which would decrease this time. 

The emerging physical picture of PuCoGas is one 
where the recoiling U nucleus generates much more dam- 
age than expected based on models of elemental Pu. This 
damage is likely dominated by near-neighbor lattice dis- 
tortions that extend into the second coordination sphere 
or beyond, possibly generating local distributions of im- 



purity phases. This damage is so severe that it encom- 
passes all the atoms in a given damage cascade. The ef- 
fective damage rate is slowed as the material anneals at 
room temperature. Photoemission results indicate that 
the 5/ electrons have both local and itinerant character 
in PuCoGa5ii£ The partially localized / electrons in the 
well-ordered material arc likely further localized in the 
damaged regions, in analogy to <5-Pu^ This view is con- 
sistent with x-ray absorption near-edge measurements on 
aged PuCoGa 5 samples^ Radiation damage therefore 
probably creates non-superconducting material both due 
to enhanced localization and strong defect scattering in 
these regions, although a proximity effect could still al- 
low some superconductivity. Annealed areas within a 
damage cascade probably would not be superconducting 
due to their limited size, unless they reach the edge of a 
cascade. 

These data support and extend the conclusion of Far- 
nan et aZ.— ^ that radiation damage occurs at a much 
faster rate than current theoretical predictions. In their 
work, the measured damage production rate is about 5 
times higher than the theoretical prediction for the un- 
relaxed defect production rate in zircons. In the present 
work, we find a production rate at least 10 times faster 
than the prediction in an intermetallic, and observe de- 
viations from a percolation model that we ascribe to an- 
nealing effects. The role of annealing should be care- 
fully considered in studies of zircons and related potential 
nuclear- waste storage materials. 



V. CONCLUSION 

Local structure data on samples of PuCoGas demon- 
strate a well-ordered local lattice structure that agrees 
with the long-range average structure obtained by diffrac- 
tion measurements. After the sample has aged long 
enough to accumulate a significant total number of a- 
decays, the local structure exhibits strong disorder, pri- 
marily through a reduction in the amplitude of the EX- 
AFS oscillations, but also in the pair-distance distribu- 
tion variances cr 2 's. This disorder affects between 20-30% 
of sample A after one year (~ 4 x 10~ 6 a-decays per 
atom), followed by a somewhat slower damage accumu- 
lation rate. Theoretical estimates that only account for 
defects predict at least an order of magnitude less dam- 
age, not including the damage-reducing factors of self- 
annealing and defect migration. These data help explain 
the fast reduction of the superconducting transition tem- 
perature both in terms of defect scattering and a simple 
percolation model. 

These results underscore the need for more local struc- 
ture studies of radiation damage in general and especially 
in PuCoGas. Only through better theoretical models and 
atomic-level probes can we understand the detailed elec- 
tronic and structural properties of damaged regions and 
how they couple to superconductivity. In particular, di- 
rect comparisons between damage cascade structures and 
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local structure measurements should be pursued. Such 
studies will have ramifications not only for understand- 
ing superconductivity in PuCoGa 5 , but also the unusual 
properties of 5-Pu, and the field of radiation damage in 
general. 
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